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The coherent-potential approximation (CPA) is extended to describe satisfactorily the motion 
of particles in a random potential which is spatially correlated and smoothly varying. In contrast 
to existing cluster-CPA methods, the present scheme preserves the simplicity of the conventional 
CPA but leads to a momentum and frequency-dependent self-energy. Its accuracy is checked by a 
comparison with the exact moments of the Green's function, and with the spectral function from 
numerical simulations. The scheme is applied to excitonic absorption spectra in different spatial 
dimensions. 

PACS numbers: 71.35.Cc, 72.80.Ng, 78.40.Pg, 78.20.Bh 



I. INTRODUCTION 

Disorder is a pertinent feature in many solid state sys- 
tems, and has been a topic of interest and intense research 
for many decades. The prototype system is a binary al- 
loy where the chemical constituents are randomly (and 
spatially uncorrelated) placed on the sites of a regular 
lattice. Adding next neighbor hopping for the electronic 
band under consideration, and leaving out any additional 
interactions, we are left with the discrete Anderson model 
which has been widely studied, mostly in view of calcu- 
lating the electronic density of states. 

A break through in the theoretical treatment was 
the invention of the coherent-potential approximation 
(CPA), initially proposed by Soven [1] and further de- 
veloped by Velicky, Kirkpatrick, and Ehrenreich [2, 3]. 
In a nutshell, the system is replaced by a perfect crystal 
(effective medium) with a complex self energy instead of 
the disorder potential, whose value is fixed by demand- 
ing that the electron at the central site is not scattered 
off the surrounding medium (on the average). The stan- 
dard (single-site) version of the CPA is characterized by a 
momentum-independent self energy, and leaves out finer 
details in the density of states. It has been extended by 
embedding not only a single site but a larger cluster into 
the effective medium. Consequently, in this embedded 
cluster CPA [4], a self-energy matrix instead of a sin- 
gle element has to be determined self-consistently. This 
can work satisfactorily for one spatial dimension [5] , but 
apart from increasing the numerical work, it is prone to 
numerical instabilities and shows sometimes an incorrect 
analytic behavior. Such problems have been reported by 
Mills and co-workers [6] when introducing the " travelling 
cluster method" for the bulk case. 

Not unexpected, a single-site approximation will fail 
completely if already the underlying disorder has a spa- 
tial correlation. This happens for an alloy with cluster- 
ing in the chemical composition. For bulk semiconductor 
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material with impurities, the finite range of the impurity 
potential itself gives rise to a spatial correlation, inde- 
pendent on the placement of the impurity atoms. Among 
the papers which have applied and refined cluster-CPA 
methods for correlated disorder, we quote Ref. [7] where 
further references on earlier work can be found. System- 
atic improvements have been achieved by Jarrell applying 
his dynamical cluster approximation (DCA) to disorder 
problems [8]. This approach has been later adopted to 
correlated disorder as well [9]. However, in order to pre- 
serve the correct analytic properties in the DCA, the mo- 
mentum dependence of the self-energy had to be coarse 
grained according to the cluster size used (see the critical 
discussion of Rowlands [10]). A cluster CPA with only 
one complex function of frequency to be determined self- 
consistently has been derived recently in Ref. [11], but 
was applied to next-neighbor correlations only. 

Our main motivation to deal with correlated disorder 
came from excitons (Coulomb-bound electron-hole pairs) 
where the averaging with the ground-state wave func- 
tion produces a correlated potential quite naturally. The 
reader who is not so much interested in exciton physics 
may skip the corresponding Sec. II and proceed imme- 
diately with the formulation of the problem in terms of 
the continuous Anderson model (Sec. III). The coherent- 
potential approximation is extended to a (smoothly) cor- 
related potential landscape in Sec. IV, introducing a mo- 
mentum and frequency-dependent self-energy. The re- 
sulting scheme preserves the simplicity of the standard 
CPA in so far as for each momentum and frequency, a 
single root searching in the complex plane is sufficient. In 
Sec. V, the exact moments of the spectral function are 
compared with those of the correlated CPA. Numerical 
results for the spectral function are presented in Sec. VI 
and compared to numerically exact simulation results. 
For the latter, the Kernel polynomial method (KPM) [12] 
is applied which rests upon an expansion of the Hamil- 
tonian into Chebyshev polynomials. In Appendix A it 
is shown that the present method guarantees the proper 
analytic behavior (Hergiotz properties) of the self-energy 
in the complex- frequency plane. Explicit expressions for 
several functions needed in the correlated CPA are listed 
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in the Appendix B while in Appendix C results for the 
exact moments are collected. 



II. APPLICATION: EXCITONS IN 
DISORDERED SEMICONDUCTORS 

Pioneering work on excitons and disorder has been 
done by Baranovskii and Efros [13] emphasizing the dis- 
tinct role of the relative motion of electron and hole 
within the exciton in contrast to the ccnter-of-mass (cm) 
motion of the exciton as a whole. Both motions aver- 
age over the underlying disorder and lead to a double- 
step "motional narrowing" which gives rise to an inho- 
mogeneous exciton line width well below the width of the 
band-edge fluctuations. In the first (relative motion) step 
one needs a proper definition of the "excitonic volume" 
which is related to the integral over the fourth power of 
the relative wave function [14, 15]. This has been refined 
later to take into account the magnetic-field dependence 
of the relative exciton motion [16]. The implementation 
of the second (cm motion) step was more complicated. 
An early attempt [17] was using a CPA-like treatment of 
electrons and holes separately but has decoupled the for- 
mation of excitons from the disorder infiuence following 
Ref. [3]. 

If the exciton is well localized on a single lattice site 
(Prenkel exciton), only the cm motion has to be con- 
sidered. This is the realm of excitations in tight-binding 
chains where disorder treatments have a long history. An 
early comparison between simulation and CPA for un- 
corrclatcd potentials has been presented by Hubcr and 
Ching [18]. The corresponding motional narrowing has 
been discussed widely based on simulations [19]. Corre- 
lated disorder of a very specific form (dimers of identical 
energy) has been treated in Ref. [20] via numerical av- 
eraging over disorder realizations. In Ref. [21], Frcnkcl 
excitons in silver-halide systems have been investigated 
using the standard CPA. Assuming still uncorrelated dis- 
order, the model has been generalized to potential distri- 
butions more complicated than simple Gaussian. 

In the case of moderate disorder (the inhomogeneous 
exciton line width being small compared with its binding 
energy), the ground-state exciton wave function for rel- 
ative motion may be treated as unaffected by disorder. 
Then, the task reduces to solve an effective one-particle 
problem for the cm motion. Still, the relative motion 
was reducing the fluctuations of the effective cm poten- 
tial (first step motional narrowing). However, much more 
important is that the resulting cm potential is now spa- 
tially correlated at least over distances comparable to the 
exciton Bohr radius - even if the underlying band-edge 
fluctuations were not correlated at all. This was our main 
stimulus to treat the exciton cm motion in a correlated 
potential. 

A straightforward way to do so are numerical simu- 
lations: For a given potential, the Schrodinger equation 
is solved for eigenvalues and eigenfunctions, wherefrom 



the spectral function can be extracted easily. Note that 
at zero cm momentum, the spectral function agrees with 
the absorption line shape of the inhomogeneously broad- 
ened exciton. A lot more effective is to solve the time- 
dependent Schrodinger equation for the Green's function 
which gives upon time Fourier transformation the spec- 
tral function directly [18, 22]. Mandatory, however, is the 
subsequent average over many disorder realizations in or- 
der to get a reasonably smooth line shape. This can be 
an expensive task limited by both computing time and 
memory resources. Therefore, it is compelling to look 
for a scheme such as CPA which avoids large scale com- 
putations but the CPA needs to be adapted to spatial 
correlation as done in the present paper. 

For excitons in a mixed semiconductor crystal, the ba- 
sic random quantity U (r) is the local energy gap due to 
A and B atoms placed uncorrelated on the lattice sites 
(see Eq. (6)). In this case, the smoothing function W{r) 
is determined by the exciton relative wave function. As- 
suming for the Is ground-state exciton the standard ex- 
ponential expression, we have 

W{r) oc £7e-''/« (1) 

where the potential correlation length ^ is proportional to 
the exciton Bohr radius Ub ■ In semiconductor nanostruc- 
tures as quantum wells (quantum wires), an additional 
source of disorder is the fluctuation of the well width (the 
wire cross section), which leads to a variation of the con- 
finement energies. If the length scale of these fiuctuations 
is smaller than the exciton Bohr radius, it can be approx- 
imated again as a delta-correlated random process U{r). 
The smoothing function W{r) is then the appropriate 
relative exciton wave function in lower dimensions. 

In the case of electrons moving in the potential of 
densely spaced randomly distributed impurities, W{r) is 
proportional to the (screened) potential of a single im- 
purity. In Appendix A, we give explicit expressions for 
exponential and for Gaussian type of the potential corre- 
lation, and are listing more details taken from Ref. [23]. 

A more sophisticated disorder-related effect is 
quantum-mechanical level repulsion of excitons which 
can be seen in time-resolved Rayleigh scattering [24] or 
spectral correlation spectroscopy [25] . While these effects 
are of two-particle nature we concentrate in the present 
paper on the one-particle level which is sufficient to de- 
scribe the absorption line shape of excitons. For a review 
on our work on excitons and disorder in nanostructures 
see Refs. [23, 26]. Quite recently, we were able to do sim- 
ulations for optical spectra without the separation into 
relative and cm exciton motion [27]. 



III. FORMULATION OF THE PROBLEM 

Let us consider a single particle (electron, exciton,...) 

with mass M which moves in a random potential V{r). 
The one-particle Green's function obeys the inhomoge- 
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neous Schrodinger equation 



IV. CORRELATED CPA 



Ar + V{r) - hz g{y, r', z) = -d{r - r') , (2) 



2M' 



with the complex energy hz. This problem is usually 
called Continuous Anderson model. The one-particle 
density of states is defined as 



dr JmC/(r, r, w — iO) 



(3) 



Further, we are interested in the spectral function Ak{u)) 

Ai,{co) = ^Q drdr' e'^^''-''">Img{r,r',uj-iO)j . (4) 

(n is a normalization volume). For cxcitons, r has to 
be understood as the center-of-mass coordinate, and the 
appropriate mass is M = rrie + mh. The line shape of the 
excitonic absorption (into the Is ground state) is given 
by the spectral function at k = ^ught ~ 0. 

The angular brackets in Eqs. (3 and 4) represent the av- 
erage over all realizations of the random potential ^(r). 
We take the average potential as zero of energy (that is 
the band edge of the virtual crystal). Therefore, the dis- 
order potential has zero mean value. Further, we assume 
it to be Gauss distributed with variance a, and define its 
spatial correlation via 

{V{r))=0, {V{r)V{r'))=a'g{r-r'). (5) 

The correlation function g{r) is therefore normalized to 
g(r = 0) = 1. All higher order correlations are assumed 
to factorize. Such type of potentials can be generated 

by integrating with a smoothing function H^(r) over a 
delta-correlated basic random process U{r) 

V{r) = j dr' W{r - r') U{r') , (6) 
Mr)f/(r')) = 6{r-r'). 

The correlation function is related to W via the convo- 
lution 



a'^g{r) = j dr' W{r + r') T^(r') . 



(7) 



Let W{r) decay on the scale ^. Then, the resulting ran- 
dom potential varies smoothly over distances of order ^, 

and the correlation of the potential can be quantified by 
the kinetic energy on the length scale ^ 



2M 



(8) 



Obviously, there are only two energies in the theory: The 
correlation energy Ec, and the disorder strength a. One 
of them can be taken as energy scale. Given a certain 
type of correlation (exponential or Gaussian), the spec- 
tral function (or exciton line shape) depends therefore 
only on the single parameter E^/cf. 



The basic idea of the CPA is sketched as follows [1, 2]: 
The entire space is divided into a central cluster where 
the disorder average has to be performed exactly, and 
into the remaining outer part which is described as an 
effective medium with averaged (but complex) potential, 
respectively, self-energy T,[z). For the determination of 
E(2;), the Green's fmiction of the effective medium G{z) 
is put equal to the cluster average of the full Green's func- 
tion (cluster embedded into the effective medium). In the 
conventional CPA, the central cluster is reduced to a sin- 
gle site, and the self-consistency condition for ^{z) can 
be formulated as a vanishing of the averaged t-matrix 
(Soven's equation, [1]). In a correlated potential, the 
abrupt change between the on-site potential ^(r = 0) 
and the self-energy T,{z) outside is in contrast to the 
smooth behavior of any potential realization. Choosing 
a larger cluster is expected to improve, but - apart from 
analyticity problems encountered - the abrupt change 
at the cluster boundary remains unaffected. In order to 
illustrate this point we look for the conditional probabil- 
ity to find a potential value Vi at a distance r keeping 
Vo = V{r = 0) fixed: 
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27rcri(r) 



cxp 



(1^1 - V^g{v)f 
2al{v) 



(9) 



The most probable value of Vi = V{r) is therefore given 
by 



V<yg{v)., 



(10) 



and distributed with a variance cri(r) = (j^(l — ^^(r)) 
increasing with distance towards cr (see Fig. 1 for a visu- 
alization). The main idea of the present CPA extension is 
to use a smooth interpolation between the central "most 
probable" potential shape and the effective medium out- 
side which is assumed to have a nonlocal self-energy, 

V{v) Vo ff(r) 6{r - r') + S(r - r', z){l - g{r')) . (11) 

In this way, any sharp break is avoided, and the effective 
Schrodinger equation reads 



/ 



h^ 

2m 



Ar + Vogir)-hz)g{r,r",z) 



(12) 



dr'i:{r-r',z){l-g{r'))g{r'y,z) = -5{r-r"). 



Fourier transformation into reciprocal space yields 

(ek + Sk(-z)-M^?kk'(^) + (13) 
+ (yo-Sk(^))^5kk^?kk'W = -^kk', 

with dispersion Ck = h^k^/2M. The Fourier transform of 
the potential correlation function, gj^, can be expressed 
via the smoothing function as 



Wk 



dre-'^^Wir). (14) 
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FIG. 1: (color online) Probability distribution of the Gauss 
correlated potential around a fixed potential value at a; = 
(top: Vb — 0, bottom: Vb = —2.5 a). The curve with dots 
marks the most probable value Vog{x). Other lines refer to 
probabilities 0.9, 0.7, 0.5, 0.3, 0.1 relative to the maximum 
(at each x value). 



The normalization ^(O) — 1 transforms into ^Zkffk — 1- 
Eq. (13) can be written as 



<5kk' + (l^o-Sk(z))V5kk^kk'( 



k 



(15) 

introducing the Green's function of the effective medium 

1 



(16) 



which is diagonal in momentum space. In contrast, the 
Green's function Eq. (15) is nondiagonal since we have lo- 
cated the central potential Vg at r = which has spoiled 
the translational invariance. Placing Vq at position rg 
leads simply to the Green's function Q{y—Yq, r'— Tq, z\ Vq) 



(for clarity, we have marked the dependence on Vq). 
Therefore, the embedding into the effective medium 
should contain two steps: (i) an integration over ro and 
(ii) an average over the (Gaussian) distribution of the 
central potential 



1 



2na 



(17) 



The self-consistency condition in the spirit of the CPA 
reads therefore 



G(r-r',z) 



dro /dT/o^(Vb)g(r-ro,r'-ro,z;Fo) 

(18) 



or after Fourier transformation 

(ekk(^;K,))v„ =Gk(z). 



(19) 



Note that only the diagonal elements of Q-w^'iz) enter the 
CPA condition. 

In order to avoid the solution of the matrix equation 
Eq. (15) we specify our Ansatz further by factorizing the 
momentum dependence as 



akk'(^;"^^o) ^Gk(z)Hk'(^;K)). 



(20) 



Then, Eq. (15) at k = k' can be solved immediately, 
and the self-consistency condition Eq. (19) converts into 
(i7k(-z; Vb)) v^, = 1. Explicitly, we have 
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(21) 



i-(yo-SkW)i?k(?i^;{s})/ 

with the auxiliary function 

i?k(?..; {E}) = E,k-.G.(.) = E . 

(22) 

The argument for chosing the present Ansatz Eqs. (11,20) 
was to exactly reproduce three well-known limiting cases: 
(i) White noise: For an uncorrelated potential (gk — 
const.), the auxiliary function Eq. (22) and therefore the 
self-energy do not depend on momentum. Thus, the con- 
ventional CPA is fully recovered which reads [2] 



l-(yo-S(z))^Ek n.^.Ui^) 



(23) 



Vo 



Here, a finite band width has to be implemented by re- 
stricting the sum over k to the first Brillouin zone, having 
J\f discrete momentum points. The parabolic dispersion 
Ek must be replaced by a tight-binding expression, and 
one arrives at the CPA of the discrete Anderson model 
with site disorder. 

(ii) Classical limit: If the correlation has infinite range, 
we have g{r) = 1 (.gk-k' = (^kkOi ^^"^ solution of 
Eq. (21) is simply 



hz- ek- l^kjz) 
hz-e-^- Vq 



= 1 



(24) 



Vo 
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with the result 

Gk(2) 



hz-Ck- Vo 



(25) 



Vo 



The Green's function mirrors the potential distribution 
displaced by ek as expected in this case. 

(iii) Second Born approximation: For weak coupling, 
Eq. (21) can be expanded into a geometrical series 



(26) 



Up to second order (j < 2), we need the disorder averages 
(Vo) = and {Vq) = cr^. This results in leading order to 
the self-energy 



z) 

(27) 

which is just the self-consistent second Born approxima- 
tion, also called "random coupling method" [28]. This 
agreement was the guiding principle for selecting the ap- 
propriate "switch on" of the self-energy in Eq. (11). We 
were starting with Vb5(r)(5(r — r') -I- I](r — r', z)(l — ^(r')) 
with a yet unspecified function ^(r'). Consequently, a 
second auxiliary function R appears as factor to the (ex- 
plicit) self-energy in Eq. (21) and Eq. (26), while Vq stays 
with the original R {R is defined like Eq. (22) but contains 
instead of g^) ■ In second Born quality, the self-energy 
reads now = a^R^/R (with arguments suppressed). 
The comparison with Eq. (27) dictates R = R that is 
g = q, which justifies a posteriori the choice made in 
Eq.(ll). 

An application of the self-consistent second Born ap- 
proximation to excitons in disordered semiconductors can 
be found in Ref. [29]. In diagram language, this self- 
energy consists of the simplest diagram where the aver- 
aged Green's fimction appears just once. Although be- 
ing an approximation with inferior quality (see Sec. VI), 
already here a self-energy comes out which depends con- 
tinuously on momentum (the frequency z plays the role 
of a parameter). Other approaches like Jarrell's DC A [9] 
- while much more sophisticated with respect to local 
correlations - fail to reproduce the self-consistent second 
Born approximation in the appropriate limit, coming up 
with a discontinuous momentum dependence of the self- 
energy. 

The coupling of self-energies at different k points in 
Eqs. (21 and (22) complicates the numerics since the self- 
consistency can be achieved only iteratively. Replacing 
Ek'(^) by T,k{z) in Eq. (22) reduces the problem to a 
simple root-searching task in the complex plane for the 
single quantity Sk(2;) (at given values of k and z) 



l-iVo- ^k{z))Rkinz - Sk(z)) 



1 



(28) 



with the auxiliary function now written as 
i?k(n^) = 5] -^'^-'^' 



hz 



Ck' 



(29) 



We call this diagonal coherent-potential approximation 
and show in Sec. V that the agreement with the sim- 
ulation results is not much sacrificed. With the self- 
consistently determined self-energy, the spectral function 
is obtained as 



Ak{(j) = Im- 



z=u>—iO 



hz-ck- T,k{z) 
and the density of states can be calculated via 



ANALYSIS OF MOMENTS 



(30) 



(31) 



Before presenting numerical results, it is useful to pro- 
vide a somewhat deeper insight into the quality of the 

CPA method presented here. We focus on the moments 
of the spectral function defined as 



M„(k) 



j^A^{uj){n^-e^r . (32) 



As a benchmark, we are going to check the moments in 
CPA against the exact ones. The first moments describe 
some basic features of the spectral function: Mq = 1 
gives the normalization, Mi is the average position with 
respect to the bare dispersion ek, and M2 is the over- 
all width. More subtle is the effect of the next moments, 
where M3 can be related to the asymmetry (or skewness), 
and to the importance of tails. A distinct advantage 
is that the moments can be evaluated exactly - as long as 
they are finite (see Appendix B). The clue is the spectral 
representation of the (disorder-averaged) Green's func- 
tion 



Tiz — Uu) 

which leads to the asymptotic expansion 

M„(k) 



Gk(^) J2 



(^^ - ek) 



(33) 



(34) 



for \z\ 00. The exact operator expression 

Gi,{z) = {k\ {hz - H)~^ ]k) (35) 

can be expanded in a similar way, and the moments are 
found as disorder averages of powers of the Hamiltonian 
in Eq. (2) 



Vo 



M„(k) = (k| {n - ek)" |k) . 



(36) 
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For the evaluation, it is useful to work in reciprocal space 
with 

(k|7^|k') = ek^kk' + Vk-k' , (37) 

where Vq is the potential of a single realization in Fourier 
space. The first steps give immediately 

Mo = l, Mi = (yk_k)=0, (38) 

since averages over odd powers of the potential vanish. 
In the next step 

M2 = 5^(Vic-k' Vk'-k) = <T^ (39) 

k' 

holds, using (VqV-qi) = a^Sqqigq and the normaliza- 
tion of the potential correlation. Interestingly enough, 
the first three moments are universal, i.e. they do not 
depend on the correlation energy, and have no k depen- 
dence either. In higher orders, a mixing between kinetic 
energy and potential correlation appears, which can be 
quantified by 

if„(k)=^5k-k'(ek'-ek)" . (40) 

k' 

For the third moment, we obtain 

Mg = a^Ki (41) 

which is again independent on momentum due to the 
mirror symmetry gq = g-q. In the next order 

M4(k) = a^K2{k) + Sct"* , (42) 

where the second term stems from the average over four 
potentials, (F^(r)) = 3cr^. We proceed up to the fifth 
moment 

M5(k) = a^Ksik) + lOa^Ki . (43) 

Explicit expressions for different correlation type and di- 
mensionality arc collected in Appendix C. 

At still larger n, things are getting involved rapidly. 
By dimensional arguments, i4'n(k = 0) oc and the 
moments at k = have the following structure 

ln/2] 

M„(0) = '^nj i?r'' (44) 

j=l 

which shows up already in the first moments listed above. 
Only the leading a power has a simple structure, namely 
a„ j=„/2 = 1 • 3 • • • (n — 3) • (n — 1) for even n. Taken alone, 
these terms would reconstruct just the Gaussian poten- 
tial distribution (classical limit). Our earlier work using 
the derivation of moments is summarized in Ref. [26]. 
The exact coefficients anj and thus the moments have 
been generated using symbolic manipulations. However, 
the numerical load increases exponentially and we were 



restricted to n < 20. On the other hand, the numeri- 
cal results gave clear evidence that much more moments 
(of the order of thousands) are needed for a proper de- 
scription of the spectral function. Consequently, earlier 
attempts to work with an infinite (but restricted) subset 
were not successful [30]. In particular at large Ec/a, even 
negative portions of the spectral function may appear. 

To complete the comparison, we are now evaluating 
the moments in CPA. First, the general relation between 
moments of the spectral function and of the self-energy 
is established. In analogy to Eq. (34), the self-energy 
moments S'i(k) are defined by 

1=0 (^^ - ^k) 

We have no c;onstant term in this expansion since the 
band edge of the virtual crystal was taken as zero of 
energy. Inserting Eq. (45) into Eq. (16) gives together 
with Eq. (34) 

ri-2 

Mo = l, Mi=0, M„(k) = ^5,(k)M„_i_2(k). 

(46) 

A similar reciirsion relation between moments holds for 
the electron gas [31] with Coulomb interaction instead 
of disorder. The auxiliary function Eq. (29) decays at 
least as l/£^ where E = 7i2; — ek in the following. Keeping 
terms up to E~^, we restrict the expansion Eq. (26) of the 
defining CPA equation Eq. (28) to j < 5 and perform the 
disorder average over the " central" potential Vq using 

{V,')=a\ (1/0^) =3^4 (47) 

(odd orders do vanish). After division with R we arrive 
at 

E = (cT^ + E2) R - Sa^-ER^ + 3a^R^ (48) 

which is accurate up to E~^. Therefore, thanks to 
Eq. (46), we can obtain the spectral moments up to M5 
as desired. In a next step, we write the auxiliary function 
needed in Eq. (28) as 

and expand this up to 1/ E^ 

^^-E + E^+ ^3 + ^4 • (5U) 

Now we are ready to equate powers of E~'" in Eq. (48) 
iteratively. In successive steps using Eq. (50) we obtain 

So{k) = 

5i(k) = a^K, (51) 

52 (k) = a^K2(k) + 2a^ 

53(k) = a^K3{k)+6a^Ki. 
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Note that the first two orders in Eq. (51) have been al- 
ready exploited to write the two last terms in Eq. (50) in 
compact form. The recursion Eq. (46) gives 

M2 = So = <J^ 

Ah = Si= a^Ki (52) 
M4(k) = S2ik) + S^ = a^K2{k)+3a^. 

Therefore, the correlated CPA generates the exact mo- 
ments up to the fourth order. Differences show up in the 
fifth order. Here, the diagonal CPA version produces 

M^^^'^ik) = a^Ksik) + 8a^K, (53) 

while the correct prefactor of the second term should be 
10 (see Eq. (43)). Using the more complicated original 
form Eq. (21) gives not a real improvement: Up to M4(k) 
there is no change but the mentioned numerical prefactor 
in Eq. (53) goes up to 9 only. 

The self-consistent second Born approximation 
Eq. (27) even fails to give M4(k) properly - it has only 
2(7'* instead of 3(t'* in the last line of Eq. (52). Using 
diagram language, it is exactly the crossing diagram in 
4th order which is missing here - while it is contained in 
CPA. 



VI. RESULTS AND DISCUSSION 

The quality of the present CPA method is judged using 
simulation results which can be considered as exact solu- 
tion of the continuous Anderson problem. For the sim- 
ulation, the r space is discretized on a cubic mesh with 
step size A, and for the _D-dimensional cube (side length 
A'^A), periodic boundary conditions are applied. In or- 
der to avoid discretization artifacts, the potential should 
change smoothly along A. We found ^/A = 3 as a suffi- 
cient condition. The standard discretization of the Lapla- 
cian operator in Eq. (2) maps the problem to the discrete 
Anderson model with correlated potential. The corre- 
sponding transfer energy is given by T ~ A^'^ /2M . A 
straightforward diagonalization would give all eigen- 
values and eigenfunctions but is restricted to unaccept- 
able small sizes A'^ because of memory size and com- 
putation time. Earlier, we had numerically solved the 
corresponding time-dependent Schrodinger equation and 
generated the spectral function by time Fourier transfor- 
mation following Glutsch [22]. Later on, we have imple- 
mented a numerical generation of Chebyshev moments, 
too [26] . Going further this way, we apply in the present 
work the powerful Kernel polynomial method (KPM) as 
detailed in Ref. [12]. In essence, successive moments of 
the Hamilton operator in terms of Chebyshev polyno- 
mials are generated. The relation to the standard mo- 
ments Eq. (32) is straightforward. However, the essential 
difference is that here not the exact (disorder-averaged) 
moments are generated but those of a given potential 
realization. The Chebyshev coefficients are damped ac- 
cording to the Jackson algorithm which gives in the spec- 
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FIG. 2: (color online) Optical density for exponentially cor- 
related disorder in two dimensions. CPA calculations in the 
diagonal version are shown for increasing ratios Ec/a (increas- 
ing peak height). To judge the quality, difFerences to the sim- 
ulation results are shown on the same scale at the bottom. 



trum a nearly Gauss-shaped line for each of the eigen- 
values. In this way, a smooth spectral function can be 
obtained after adding up results of a sufficient number 
of independent disorder realizations. We have carefully 
checked that all these technical parameters are chosen 
such that no influence on the final shape of the spectral 
function is seen. For a typical calculation in _D = 2, a 
square grid with N — 100 has been used, and an average 
over 3000 realizations was performed. In order to get a 
reasonably smooth spectrum without too much broad- 
ening, 500-1500 Chebyshev moments were calculated for 
the Jackson algorithm. The numerical effort was maxi- 
mal in D = 3, where 500 realizations for a box with side 
length A^ = 50 have been added up. 

We begin with a comparison of results in D = 2 which 
are relevant for quantum wells with disorder. The cor- 
relation type is taken exponential, and we concentrate 
here and in what follows on the spectral function at k = 
which gives directly the inhomogeneous broadening of the 
exciton absorption line, called optical density in the fol- 
lowing. Starting with the Gaussian potential distribution 
in the classical limit [Ed a = 0), the curves in Fig. 2 are 
getting more narrow and asymmetric for increasing val- 
ues of Eel (7. This represents the motional narrowing on 
the cm level. The correlated CPA deviates only slightly 
from the numerically exact results, as visualized by the 
bottom curves. 

In Fig. 3, we compare the case Ecjcr — 2 for differ- 
ent levels of approximations. The self-consistent second 
Born approximation (dashed-dotted) deviates markedly. 
Only its high-energy asymptotics is reasonable since it is 
dominated by perturbation theory. In particular, it fails 
completely in the low-energy tail where a sharp (square- 
root) cutoff is produced. This is an inherent feature of 
any diagrammatic expansion which leads to a geometri- 
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FIG. 3: (color online) Optical density for exponentially cor- 
related disorder in two dimensions at Ec/o — 2. Different 
levels of approximation are compared with the simulation re- 
sults (see text). 



cal series in terms of the interaction (disorder). The CPA 
does here pretty well since the final summation over the 
local potential fluctuations brings in a true random fea- 
ture. The full version (dotted fine, Eq. (21)) gives only 
a slight improvement compared to the diagonal approx- 
imation (full line, Eq. (28)). However, the numerics for 
the full problem is much more involved since at a given 
energy, the complete momentum-dependent self-energy 
has to be brought to convergency. We succeeded only by 
using acceleration, respectively, slowing down in the re- 
cursive determination. On the other hand, the diagonal 
version needs only a single zero search in the complex 
plane. Therefore, we show in all the other figures exclu- 
sively results from the diagonal version. Then, for the 
optical density, only I]k=o(z) has to be determined self- 
consistently. 

It is pleasing to see how the quality of the differ- 
ent approximations (self-consistent second Born, diag- 
onal CPA, and full CPA) goes in parallel to the num- 
ber of moments which are exactly reproduced. Taking 
all other parameters unchanged, a correlation of Gauss 
shape leads to narrower spectral functions compared to 
the exponential type (Fig. 4(a)). The low-energy tail is 
shown on a semi-logarithmic plot in Fig. 4(b). The sim- 
ulation is getting noisy there since states deep in the tail 
are rare events but compares well with the asymptoti- 
cally strict result of the optimum fluctuation theory [32], 
extended here to finite spatial correlation [26] . While the 
correlated CPA follows initially rather close, deep in the 
tail the spectral function is somewhat overestimated. 

When going from two to three dimensions (Fig. 5), 
the line width at a given value of Ec/a is reduced (cp. 
Fig. 4(a)). This can be related to the wave-function ex- 
tension as - loosely speaking - an increased localization 
length. As well known, it is harder to localize the wave 
function in _D = 3 compared with D = 2, assuming 




-3 -2 -1 
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FIG. 4: (color online) Optical density for Gauss correlated 
disorder in two dimensions. In addition to the linear display 
in (a), the logarithmic display in (b) emphasizes the tail re- 
gion. The optimum fluctuation result (chained curve) has 
been adjusted vertically. 
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FIG. 5: (color online) Optical density for exponentially cor- 
related disorder in three dimensions. 
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FIG. 6: (color online) Optical density for Gauss correlated 
disorder in one dimension. 

the same disorder strength and correlation type. Since 
we are here interested in states around the lower band 
edge (these have the dominant contribution to the k = 
spectral function), the localization edge of the Anderson 
model in £) = 3 is of no relevance here. 

To complete the analysis, we show in Fig. 6 results for 
the one-dimensional problem with Gauss correlated dis- 
order. A final plot (Fig. 7) deals with the uncorrelated 
case in Z? = 2 which is the realm of the standard single- 
site CPA in the discrete Anderson model. The dispersion 
to be used here 

2 

Ek = 2r^[l-cos(fcjA)] (54) 

refers to a simple cubic lattice {Af = il/A^). The de- 
viations from the (numerically exact) simulation results 
are definitely smaller than in the correlated cases (Figs. 2 
and 4). This is in complete accordance with the mo- 
ment analysis, since the standard tight-binding CPA for 
uncorrelated disorder preserves the exact moments even 
up to Mj(k) (we quote the first orders in Appendix C, 
Eq. (C7)). Still, the CPA extension towards correlated 
potentials is of reasonable quality, and provides therefore 
a new tool for studying disorder problems in solid state 
physics. 

The present method could be used for the discrete case 
with spatial correlation as well. The potential generation 
Eq. (6) has to be discretized as 

y(R) = ^ ^^(R - R') C/(R') , (55) 

R' 

(C/(R) C/(R')) - W 

where R denote the lattice points. All Fourier trans- 
forms are restricted to the first Brillouin zone, and in 
Eqs.(22) and (29) Ck has to be taken as tight-binding 
dispersion again. However, an application to a binary 



FIG. 7: (color online) Optical density (or zero momentum 
spectral function) in a two-dimensional tight-binding model. 
The standard CPA result for uncorrelated site disorder is 
shown and compared with the simulation (here using 2048 
Chebychev moments). 



alloy with spatial correlation (clustering) of the chemical 
species is not possible: In our method, the (local) po- 
tential must be Gauss distributed due to the averaging 
process Eq. (55), while the relevant correlated potential 
for the alloy should be still a binary quantity. 
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APPENDIX A: ANALYTIC PROPERTIES OF 
THE SELF-ENERGY 

The proper Green's function of the averaged medium 
G]t:{z) - considered as a function of the complex- 
frequency argument z - should obey the following ana- 
lytic properties (called Herglotz properties, HP [6]): ft is 
analytic everywhere outside the real axis, obeys the mir- 
ror symmetry Gk(z*) = 6^(2;), and its imaginary part is 
non negative in the lower half plane. Using physics lan- 
guage, HP guarantee that the Green's function refers to 
a causal system, and the spectral function (the jump on 
the cut along the real axis, see Eq. (30)) is non negative, 
as appropriate for the probability amplitude for adding 
a particle with momentum k to the system. Due to the 
simple relation to the self-energy Eq. (16), Sk(^) should 
have HP as well. An alternative formulation is (i) the 
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validity of the spectral representation 
^k[z) - 



(Al) 



with (ii) a non negative spectral density Sk(w) > 0. 
Strictly speaking, a function with HP could have - in 
addition to the integral in Eq. (Al) - a constant and a 
term linear in z. However, in the present case such terms 
are absent since the self-energy has the virtual crystal 
as reference. Formally, they are missing in the asymp- 
totic expansion Eq. (45) as well. In the following we will 
show that the present extension of the CPA to correlated 
disorder generates a self-energy with Herglotz properties. 

The mirror symmetry T,^(z) = Ek(z*) is obvious since 
- apart from I^uiz) and z - only real functions (ek, 5k, 
P(Vo) enter the CPA equations. To be definite we place 
now z into the lower half plane (Jm z < 0) and search 
for a self-energy with ImEk(2) > 0. Then, it follows at 
once that the auxiliary function Rk{hz, {Y,}) = R de- 
fined in Eq. (22) has a positive imaginary part, i?2 > 
(note that (?k-k' is non-negative, Eq. (14)). The self con- 
sistency condition Eq. (21) is written explicitly as 



l/R + T.~ Vb 



= R 



(A2) 



where E = Sk(2^)- The integral looks like a standard 
spectral representation, but the imaginary part of the 
denominator D = \/R + Yj — Vq which is I?2 = S2 — 
R2l\R\'^ may cross zero at some curve in the lower half 
plane. This would give rise to a non-analytic self-energy 
there but can indeed never occur: The imaginary part of 
Eq. (A2) reads 



-D2 / dVr 



PjVo) 
\D\^ 



(A3) 



with a strictly positive value of the integral since P{Vq) is 
a positive probability distribution. Therefore, a solution 
can never have D2 — since i?2 > as shown above. 
This completes the proof that the self-consistent self- 
energy has HP. The diagonal simplification of Eq. (28) 
does not change any step of the proof. In this case, the 
even stronger assertion D2 < hz2 < holds since the 
inequality Im(l/i?k(^) — ^z) < can be proven quite 
generally {Hz = hz — T,k{z)). 

In order to see the desired behavior in the numerics as 
well, we have been searching in the full complex z plane 
for self-consistent solutions of the self-energy. Indeed, the 
self-energy was found to vary smoothly as a function of z 
(crossing non- analytic points/curves would show up as a 
discontinuity). Still, there could be more than one solu- 
tion. Assuming at large \z\ a, vanishing self-energy as ini- 
tial guess in the root search, we were at least starting with 
the proper solution, and have continued with the previ- 
ous solution as start for the next z position. A branching 
of this solution into two is not possible - this would signal 
a non-analytic point. However, there could be an acci- 
dental degeneracy with a second solution. Although this 




FIG. 8: (color online) CPA self-energy for exponential corre- 
lation {Ec/(J = 1 in two dimensions) at zero momentum (top) 
and at Ek = o" (bottom). The dots give the real part calcu- 
lated via Eq. (A4) . The exact self-energy for ek = from the 
simulation is shown dashed. 



never occurred in our calculations, we could imagine how 
to select numerically the proper continuation (demand- 
ing the absence of breaks in slope) . We were not able to 
give a general proof for the uniqueness of the solution, 
as done by Mills and co-workers in Ref. [6] for the substi- 
tutional alloy with uncorrelated disorder, and in Ref. [33] 
for a chain with randomly placed delta scatterers. A 
more practical proof for Sk(z) having HP is to check that 
the generated function obeys the spectral representation 
Eq. (Al). To do so it is sufficient to run the calculation 
for z = Lu ~ iO only, since Sk(i^) = -rmSk(^ — *0). The 
result in Fig. 8 shows a complete agreement between the 
direct result and the spectral form 



/oo J f -| 
— Sk{u')V 
-00 UJ-UJ' 



(A4) 



where V stands for the principal value. For compari- 
son, we have added in the figure exact self-energy results 
which follow from the complex Green's function gener- 
ated with our simulation technique using Eq. (16). 



APPENDIX B: EXPLICIT EXPRESSIONS FOR 
DIFFERENT CORRELATION TYPES 



The probability distribution of the local potential value 
Vq is of Gauss type, Eq. (17). Therefore, the integral in 
Eq. (A2) can be expressed using the complex error func- 
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tion 



resulting in 



iu{z) = 



z-t 



(Bl) 



(B2) 



For the diagonal version, R = R]s_{hz — Sk(-z)) is un- 
derstood. The specification of the potential correlation 
function depends on the system under investigation. If 
the potential is short-range correlated, wc approach the 
limit of a "white noise potential" (gk = 1/^)- This is the 
proper situation for the traditional CPA where the self- 
energy is assumed to be site diagonal (i.e. independent 
on momentum k). 

For the exciton case, it is important to note that elec- 
tron and hole are scanning the potential landscape on a 
different scale, depending on their masses m,. and m^j. 
For quantum wells, local fluctuations of the well width 
Lz lead to local shifts of the band edges Ea{Lz), and the 
smoothing function has the following form [23] 



W{r) = hC 



a=e,h 



dEg 
dL, 



(B3) 



with mass factors -qe = M/m,h and rjh = M /m,,. The pa- 
rameters h (C) characterize typical height (lateral exten- 
sion) of the well width fluctuations (island size) . The ex- 
citon wave function can be taken hydrogen like, (/'is(r) cx 
exp(— r/as), where as is he appropriate Bohr radius of 
the quantum well exciton. For equal electron and hole 
mass (?7e = r]h = 2), Eq. (B3) reduces to a single expo- 
nential dependence 



W{r) = 



2 a 



-r/i 



where (, = as/4 and 



(B4) 



(B5) 



Performing the two-dimensional integrations we evaluate 
Eq. (14) with the result 



9k = 



1 



(B6) 



Roihz) 



dx 

wJo {nz-EeX^){l + x^)^ 
s + 3- 21 
E,(\ + s)2 • 



(B9) 



Finally, we apply Eq. (1) to the bulk mixed crystal and 
get 



5k = 



647r^3 Y 



n (1 + {ki)2y 



which is followed by 



Roihz) = — [ 

Jo 



oo 2 



x'^dx 1 



(BIO) 



(Bll) 



'o hz-E^x"^ (1-Fa;2)4 

_ + + 15s - 5 + l&sf^ 

~ Ec{l + sY ■ 

For a Gauss-type correlation with characteristic length I, 
we assume W{r) cx exp(— r^/2£^) and obtain 



(B12) 



where D = 1,2,3 is the spatial dimension. The correla- 
tion energy is now defined as E^ = Ti^ /{2Mf'), and we 
have to evaluate 



Roifiz) 



1 

ErT^D/2 



s — x-^ 



(B13) 



which gives 



n {1 + [k^YY • 

The auxiliary function R\^{?iz) Eq. (29) at k = reads 

dx 2 s"^ +4:3 + 3 + 2log{-s) 



D = l : Ro{nz) = i^^w{-V^), 

D = 2 : Roihz) = -^Ei{-s) , (B14) 

-C'c 

D = 3 : Roihz) = ^ (^^^^^(;(-^/i) - l) . 

In addition to the complex error function Eq. (Bl), the 
exponential integral E^iz) = dte~*/t enters. 



APPENDIX C: RESULT FOR THE EXACT 
MOMENTS 

We list here explicit values for the first exact moments 
Eq. (36). For Gauss correlation, we take advantage of the 
closed form 



f 

Roihz) = / 
Jo 



hz - EcX il + a;)3 



E4l + s)^ 

(B7) 

where s = hz/ E^ is the dimensionless complex energy. 

The corresponding results for exponential correlation 
in quantum wires (one dimension) are 



K„(k = 0) = £, 

in D dimensions, and obtain 
D 



„rCD/2 + n) 
ViD/2) 



(CI) 



(C2) 



5k = 



4^ 1 



(1 + (fcO')' ' 



(B8) 



M4(k) 

M5(k) 



a'Ei 



iD + 2)D 



2iklf 



+ 3a^ 



iD + 4)D 



+ Qiklf 



+ 5Da'^Ec . 
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In the exponential case 
Ms 



(C3) 



holds independent on spatial dimension. The last finite 
moment is here 



M4(k) = a^E^^ 



4 9' 
3«' 



+ 3(7^ 



(C4) 



for D = 3. This unexpected termination of the moment 
expansion can be understood quite easily: The spectral 
function Eq. (30) decays at large positive energies as 



^2 



(C5) 



which follows from plain perturbation theory using 
Eq. (27) in leading order. For the Gauss correlation 
Eq. (B12), multiplication with any power of hw = Ck 
leads to a convergent frequency (better momentum) in- 
tegral. For the exponential case, however, we have 



gk (x k (2^+2)^ and the integrand of the nth moment 
behaves as /c^^""-^"^). Therefore, the moments are finite 
up to n = 3 (£) = 1, 2) or n = 4 (£) = 3). 

For the sake of completeness, wc quote the moments 
in the uncorrclatcd (tight-binding) case as well. Mq = 1, 
Ml = 0, and M2 = cr"^ hold as before. We write the 
(simple cubic) dispersion as 

1 ^ 

ek = 2TD (1 - C(k)) , C(k) = cos{kjA) (C6) 



and obtain instead of Eqs. (41)-(43) 

M3(k) = 2TDa'^C{li), 

M4(k) = i2TDfa^ (c\\c) + ^) + , 



(C7) 



M5(k) = C(k) I (2^)3^2 \^C^(k) + ^]+ IQTDa^ 
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